Nonlinear Schrodinger Equation with Random Gaussian Input: 
Distribution of Inverse Scattering Data and Eigenvalues 
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We calculate the Lyapunov exponent for the non-Hermitian Zakharov-Shabat eigenvalue problem 
corresponding to the attractive non-linear Schrodinger equation with a Gaussian random pulse as 
initial value function. Using an extension of the Thouless formula to non-Hermitian random oper- 
ators, we calculate the corresponding average density of states. We also calculate the distribution 
of a set of scattering data of the Zakharov-Shabat operator that determine the asymptotics of the 
eigenfunctions. We analyze two cases, one with circularly symmetric complex Gaussian pulses and 
the other with real Gaussian pulses. We discuss the implications in the context of information 
transmission through non-linear optical fibers. 



I. INTRODUCTION 

One defining development in telecommunications tech- 
nology during the last two decades has been the 
widespread use of optical fibers for transmitting enor- 
mous quantities of data across large - even transoceanic - 
distances. For such increasingly large distances, the non- 
linearities in the fiber cannot be neglected, as they tend 
to distort transmitted pulses. Consequently, the detec- 
tion of traditionally modulated signals becomes problem- 
atic. For fibers with negative group velocity dispersion 
(GVD) it is possible to compensate these effects by cre- 
ating stable solitonic pulses [ll, 0] ■ As a first approxima- 
tion, these solitary waves are solutions of the non-linear 
Schrodinger equation (NLSE), the effective equation de- 
scribing propagation of light in the frame comoving with 
the mean group velocity 0]. In normalized units the 
NLSE is expressed as 



.du d'^u 



+ — + 2\u\'u = 



dx dt^ 



(1) 



where u(t, x) is the (complex) envelope of the electric 
field, carrying the transmitted information signal along 
the fiber [i^. Traditional analyses of this equation fo- 
cus on single and dilute solitonic propagation Q. How- 
ever, to address the ultimate information capacity limits 
through the fiber using solitonic pulses, one needs to ex- 
plicitly consider dense soliton systems, where the soliton 
interactions can no longer be treated as small. 

The problem of determining the spatial evolution of 
an incoming pulse u{t) = u{t,0) is solved via the inverse 
scattering transform (1ST), where u{t) enters as the "po- 
tential" in a linear eigenvalue problem. For the NLSE 
this is the Zakharov-Shabat (ZS) eigenvalue problem Q, 
comprising of a 2 x 2 system of coupled first order differ- 
ential equations. 
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*,(t) = ?7(i)*,(t) =^*,(t), (2) 



where ^z{t) ~ [V'i(0 "02 (i)] , and appropriate asymp- 
totic conditions on the eigenstatcs. given in the next sec- 
tion. 



In this paper we analyze the distribution of the scat- 
tering data, i.e. the average density of states (DOS) of U 
and the average distribution of a set of complex numbers 
{bz} that determine the asymptotics of the cigcnstates, 
when u{t) is drawn from a zero- mean, (5-correlated Gaus- 
sian distribution, describing the distribution of transmit- 
ted codewords. Gaussian input signals are often used 
in information theory, and in linear transmission prob- 
lems they often reach the Shannon capacity Q • In addi- 
tion, when the characteristic signal amplitude uq is much 
smaller than its bandwidth (but with D = UqT ar- 
bitrary), it is reasonable to approximate Q the input 
distribution with a (5-correlatcd Gaussian for eigenvalues 
z small in the scale of r^^. 

The non-hermiticity of U causes the eigenvalues to 
spread over the complex plane. This generally makes 
the exact calculation of the DOS more difficult. Several 
powerful methods have been developed for calculating the 
statistical properties of non-Hermitian operators, which 
appear in the modelling of diverse physical processes (see 
e.g. U^MMM,M,EA, 1^ M IT, 18, 19, 20, 21.]). In 
most cases the random matrices are treated in a mean- 
field sense and are thus considered full random matri- 
ces. However, to our knowledge there are only few non- 
Hermitian operators with diagonal randomness for which 
the exact density of states has been calculated in closed 
form p^. [23l. [2^ . In our case, we first calculate the Lya- 
punov exponent in closed form taking advantage of its 
self- averaging properties. Combi ning this with a gener- 
alization of the Thouless formula [23 for non-Hermitian 
operators [2^, that relates the Lyapunov exponent with 
the DOS, we arrive at an explicit expression for the lat- 
ter. Since the Lyapunov exponent is simply related to 
the localization length, it also provides information for 
the eigcnfuctions of U. 

In addition to the DOS we calculate the limiting distri- 
bution of the scattering data coefficients {b^}, which de- 
pends strongly on the input distribution of u(t): For cir- 
cularly complex u{t) the distribution of In bz approaches a 
Gaussian distribution albeit with singular variance grow- 
ing as TlnT, while for real u{t) the distribution is highly 
singular, approaching a Cauchy distribution. 

It should be noted that the Hcrmitian "counterpart" 
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of this operator, 



Uh 



idt u*{t) 
u[t) -idt 



(3) 



arises in the 1ST for positive GVD, and also as a special 
case of the fluctuating gap model of disordered Peierls 
chains (see [l^] and references therein). Its DOS and 
localization length have a long history of analysis [1, [28l . 

in,!!!. 

The spectrum of U , together with the asymptotic be- 
havior of the corresponding eigenstates ^ z, which as we 
shall see is determined by have the same informa- 
tion content as the input signal u{t). This is because 
inverse scattering transform mapping between the scat- 
tering data of all eigenstates and u{t) is one-to-one @ 
[1^. However, while the spatial evolution of u{t,x > 0) 
and the eigenstates ^{t,x > 0) is quite complicated, the 
eigenvalues z oi U remain constant as the signal prop- 
agates down the fiber, and the corresponding scattering 
data vary in a trivial manner In fact, they can both 
be seen as playing the role of "action" variables changing 
adiabatically in the presence of non-integrable perturba- 
tions. Therefore, the problem of light propagation in the 
fiber becomes easier to analyze in terms of the scattering 
data of the Zakharov-Shabat eigenproblem, especially in 
the presence of perturbations to ([T]) , such as noise due to 
amplification or phase conjugation, which will ultimately 
determine the optical fiber capacity [Sll, IH, [H, HI] . As 
a result, the description of the scattering data as a func- 
tion of the input signal u{t) may provide a framework for 
understanding the ultimate limits of information transfer 
through optical fibers. 



II. LYAPUNOV EXPONENT AND DOS 

We will now describe the basic steps to calculate the 
Lyapunov exponent of t/ in ([2) which will then lead to 
the average DOS. To proceed, we start by introducing 
the ZS eigenvalue problem. Traditionally, this is defined 
as a scattering problem of the operator [/ in ([2]), in the 
presence of the potential u{t) = u{t,x = 0), which de- 
cays sufficiently fast for t ±oo. In this context the 
scattering states are set up with the following asymptotic 
conditions outside the range of the potential: 



as t 



as t 
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with the a, &'s being the transmission and reflection co- 
efficients respectively. By taking into account the sym- 
metry of the problem under complex conjugation it is 
possible to show that a{z*) = a*{z) and b{z*) = -b*{z), 
where the star (*) denotes the complex conjugate. 

When the above solutions correspond to a localized 
eigenfunction with eigenvalue z, the transmission coeffi- 
cient a{z) has to vanish at that z, making the two sets of 
solutions directly proportional: 



*,(0 = -6**.(<) 



(6) 



where $z and $z are the admissible exponentially de- 
caying eigenfunctions for Im{z) > and Im{z) < 0, 
respectively. Note that inside the region where u{t) is fi- 
nite, they should decay with a Lyapunov exponent k{z), 
rather than with |/m(z)| as in The proportionality 
constants bz in ^ are not simply related to the func- 
tions b{z) evaluated at the eigenvalue z [i^. It is clear 
from above that delocalized states can only exist when 
Im{z) = 0. 

The proportionality factors bz and their corresponding 
eigenvalues z are very important quantities in the the- 
ory of the Inverse Scattering Transforms: Together with 
the continuum delocalized states characterized by b{z), 
they can completely reconstruct the original u(t). There- 
fore, in the context of information theory, they carry the 
same information content. In physical terms, the local- 
ized eigenstates of the Zakharov-Shabat problem corre- 
spond (through the 1ST) to the solitonic excitations in 
the fiber, while the continuous spectrum for Im{z) = 
gives the radiation modes, which spread out and decrease 
in amplitude as the signal propagates down the optical 
channel. We will focus on the localized states, since in 
the limit T — > oo they correspond to the dominant part 
of the solution. 

Our computation of the DOS of the problem is based 
on the calculation of the Lyapunov exponent ^(^,77), 
which then yields the density of states through the gen- 
eralized Thoulcss formula (derived in Appendix A): 



The (upper) Lyapunov exponent is defined by: 



1 

2^ 



1 



(4) which can also be written as: 



(7) 



(8) 



(9) 



For concreteness, we express the eigenvalue 2; as z = 
^ + irj. The two sets of solutions in ^ are linearly related 
through the S'-matrix: 



(5) 



Since the system is self-averaging (the evolution of ^pl , ip2 
along t is a Markov process) , we can exchange the average 
over t in ([9]) with an average over the Gaussian ensemble: 



K ^ — lim 

2 t — >oo 



|ln(|^ip + |^.n 



(10) 
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This is our starting point for calculating k. From ([2|) 
we find: 



at(|Vi|' + l^2p) =2,7(|V'i 



(11) 

„u)(t)+»0(t) 



Defining the complex variable /(t) = -^Jt^jj = e" 
with w G (—00, 00) and (f> G [0, 27r), wc can rewrite 
as: 



atln(|^i|2 + |^2p) =27ytanhw;. 



(12) 



We are interested in the long-time behavior of (|12p . For 
a given u{t), w(t) undergoes constant change at any t, 
but the probability distribution of its values in the Gaus- 
sian ensemble will tend to a stationary distribution for 
large t. To see this, we must derive the Fokker-Planck 
equation for the joint probability distribution P{w, 0; t). 
This is straightforward for (5-correlatcd Gaussian poten- 
tials, since in this case w{t) and (j){t) become Markov 
processes [ssj . 



A. Circular complex Gaussian potential u{t) 

We start by calculating the density of states (DOS) 
p{^,r]) and locahzation length l{^,r]) when u{t) is cir- 
cularly symmetric, i.e. u{t) = {ui{t) + iu2{t)) with 

ui,U2 real, {ui{t)) = 0, and {ui{t)uj(t')) = DSijS{t — 
t') , i,i ~ 1, 2. In this case, the evolution of w an " 
described by the set of stochastic equations: 

dtw — 2r/ -|- 2 coshw (sin (^ui cos 0U2) 
dtii = — 2^ — 2sinhu' (cos(/)Ui — sin(/)U2) • 

The Fokker-Planck equation derived from these (i 
Stratonovich picture) is: 

dtP = 9„ {{\-v^){-2T]^Dd^)\ P 



d,, P 



1 



where v = tanhw. A simplification can be obtain 
integrating over (j). Because the right hand side o 
depends only on v, we only need P{v) = J^^ d(j)P 
to calculate the average. Integrating over (j) and 
the periodicity of P in this variable we find the F( 
Planck equation for P: 



□ 



dtP 



[(i-i.2)(-27, + z?a„)] p. 



(15) 



Setting the left-hand side to zero we find the stationary 
solution to which the system relaxes for large t: 



P{v) 



•qe o 



Dsinhf 



(16) 



This is also a stationary solution of the full Fokker-Planck 
equation (|15p . implying that asymptotically becomes 



uniformly distributed. We can now calculate the Lya- 
punov exponent from equations PTT]) . (fT^ : 



K = r/ / dvvP{v) 



D 
"2" 



27] 
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coth 



(17) 



Note that for large \rj\, k sa \r]\ independently of D: this 
is expected since in this limit the potential decouples the 
left (-01) from the right moving {1P2) wavefunctions. A 
simple application of the Thouless formula ([7]), gives the 
exact density of eigenstates for the system: 



2 !?coth(ga)-i 
ttD sinh^ (f ) ■ 



(18) 



The independence of p from ^ is not surprising: the den- 
sity of states of the Hermitian (diagonal) part of ^ is 
independent of f. Therefore, in the so-called mean-field 
approximation [l^ . [36| the extension in the imaginary 
axis will be ^-independent. It should be noted however 
that that mean-field approach would have given a con- 
stant DOS within a zone around 77 = 0, rather than (jlSp . 
A comparison of this expression with the result of numer- 
ical simulations can be seen in Fig. [TJ Again note that 
for large rj the density of states vanishes: in this limit 
there is an exponentially small probability for finding a 
potential deep enough to create a bound state. 




FIG. 1: Theoretical curve (solid line) and results of numerical 
simulations for the profile of the DOS vs. 77. We have used 
the modified Ablowitz-Ladik diagonalization scheme [s^] to 
ensure that the determinant of the discretized transfer matrix 
has unit value. The value of D is 1, the size of the system 
is T = 135 and the step size is 0.075. The disturbance near 
77 = is a finite-size effect. The localization length grows as 
I ~ 3D/2r]^ near 77 = and so numerical results differ from 



our T ^ 00 formula for 



< 



4 



The localization length rj) is the inverse of the Lya- 
punov exponent, I = k~^. To see this, we note that the 
Wronskian of two independent solutions of ([2]) is con- 
stant, therefore if for a given z it has an solution increas- 
ing exponentially as exp(Ki), its other solution has to be 
exponentially decreasing as exp(— Kt). Thus a square in- 
tegrable solution necessarily decays with length-scale k^^ 
inside the support of u{t). From pT|) one can see that 
states become increasingly delocalized as the eigenvalues 
approach the real axis on the complex z plane: I diverges 
as I ~ ^ near the real axis. The localization length also 
determines the stability of the corresponding eigenvalue 
to the presence of a finite time window of the pulse T. 
Specifically, the typical lifetime of a state with eigenvalue 
z will scale as ~ e'^i^^'^/'^ Q. Indeed we see this in Fig.[Tl 
where for states with localization length comparable to 
the system length T, i.e. close to i] = 0, the calculated 
DOS is no longer valid. To capture the behavior of the 
DOS in this region, a zero-dimensional analysis similar 
to [23, SO, Sa, Syl is needed. 



B. Real initial pulse u{t) 

We can also analyze the opposite case when u{t) is real, 
Gaussian with {u{t)) = 0, and {u{t)u{t')) ^ DS{t - t') . 
In this case, the evolution of w and (f) is described by 
(fT3| by setting U2{t) = Im{u{t)) = 0. The corresponding 
steady state solution of the Fokkcr-Planck equation can 
be derived from: 







[(1 - t.2)(sin2 (l)Dd^, - 277)] P + 2id^P (19) 



sin 24>d^ 



1 



2vd„ 



P 



2v^ cos 
1 - i;2 



-d'iP 



For large |^|, P is independent of (p to leading order in ^. 
Therefore the large-^ expansion is essentially identical 
to a Fourier expansion. Integrating (|19p over (j) gives 
(fT5)) . Thus P is to leading order identical to that of 
the circularly symmetric complex u. After some algebra 
one can derive the next-leading order result. To order 
0(1/^^) the correction to the Lyapunov exponent is 



2if 



(20) 



resulting in the following correction to the DOS expres- 
sion of (UHl) 



Sp 



jj2 Srycoth 



D sinh^ ( 



4if_ 
D 



sinh^ (^) 



(21) 



In the opposite limit of small ^, we expect the distribution 
in (f> to be peaked. Indeed for ^ = 0, (fTg]) has a solution 
that is proportional to S{cos(j)). This results in 



Piv) 



(22) 



with corresponding Lyapunov exponent 

Vhiv/D) 



K{ri) 



loiv/D) 



(23) 



where /o,i are modified Bessel functions of the first kind. 
Wc see that compared to , is more singular when 
\v\ w 1, i.e. for large w. 



III. DISTRIBUTION OF 

The complex numbers bz, that determine the asymp- 
totics of the bound states of U, can be expressed in terms 
of the limiting behavior of the eigenfunctions. Specifi- 
cally, for 77 > we have from (|4|6p 



Defining ^z{t) = ^z(—t) we can write: 

Mm bit), bit)^t^, 



(25) 



where for convenience we have dropped the subscript z. 
The time evolution of In b is found from (O , 



= i(uf + uf 



<91n6 
dt 



with f{t) = 1^, fit) = u{t) ^ u{-t), and 



(26) 



^ = -2izf + iu*-iup 



|[ = ~2izf -iu + iu*p. 



(27) 



A. Circular complex Gaussian u{t) 

For a circularly symmetric complex Gaussian u, the so- 
lution of the Fokker-Planck equation derived from (|26l27p 
relaxes for large times towards a stationary solution 
where In 6 is uniformly distributed [4^ . while / and /, 



are 



expressed in polar form, i.e. / = e™^**", / = e 
distributed independently according to the steady state 
solution (fT6)) . Because of the infinite range of the real 
part of In b however, this stationary solution is ill-defined. 
A better approach is to discretize the size T of the pulse 
into steps of size r, equal to the inverse bandwidth of the 
input signal. Equation (|26p then reads 



(28) 



T/2r 

In 5 = ir^ [uifi + u*f^ 



The variables w^, vn are i.i.d. Gaussian random variables, 
distributed according to 



Pu{u) 



1 



2'KUo 



(29) 
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where Uq = D/2t and u stands for the real or imaginary 
part of either variable. For large enough T, the sum in 
([28|) will be dominated by the domain where the distri- 
butions of fi and fi have reached their steady state. In 
this domain, we find that the real and imaginary parts 
of the products Xi = Uifi and Xi = u*fi have zero mean 
and the tails of their distributions fall off as the inverse 
third power of the argument. More precisely. 



1 



'i?sinh(f ) 



(30) 



where x stands for the real and imaginary parts of Xj , Xj . 
The general theory for sums of random variables [40, |4lj, 
I42I then tells us that for large T/r the distribution of 
In |6| will be Gaussian, with zero mean and variance 



^Tln^. 
sinh(§) 2t 



(31) 



The imaginary part of In 6 is an angle and so, although it 
follows the same distribution as the real part, will due to 
periodicity become imiformly distributed in [0, 27r). As 
seen in (|6]) for 77 < the corresponding hz is replaced by 
—h* [4^. Thus their distribution will be the same as that 
of the &'s, with -q replaced by — r/ in ([51]) . 



B. Real Gaussian u{i) 

In this case, the Fokker-Planck equation for In 6, de- 
rived from ()26|27p after setting U2 and 112 to zero, again 
predicts that its distribution becomes uniform as the du- 
ration T of the pulse grows to infinity. Equation (|26p can 
be written as: 



dt 



= — ue sm ( 



ue cos ( 



f ue™ cos 



(32) 



As seen above, an exact solution to p9|) is not available, 
but we can still obtain the first terms of an expansion of 
the stationary probability distribution in powers of e"™: 



P{w,(j>) « aie~'"(5(cos0) -I- 



^26 



(33) 



with an identical expansion for the distribution of li and 
(p. The constants ai,a2 depend on £,/D and ij/D, but, 
being related to the normalization, they cannot be de- 
termined without a knowledge of the full solution. We 
can thus only partially specify the manner in which the 
real and imaginary parts of In b approach uniformity as 
T grows. 

As in the complex case, the real part, ln|6|, will be 
a sum of independent variables Xi = Uie^^ sin(f>i. How- 
ever, in this case, due to the more singular behavior of 
P{w, (p) for large w, the tails of Xi will be longer, faUing 
off as 1/lxiP for large T. As a result, the distribution 
of ln|6| will asymptotically follow a Cauchy distribution 



scaling like T/r [4l|, [4^. Its statistical median will be 
zero by symmetry, coming from the even parity of the 
Gaussian distribution of u, ii. The phase of b docs not 
get contributions from the first term in (j33p because of 
the delta function in this term. For large T, the second 
term in the expansion dominates, making it uniform over 
[0, 27r), in the same manner we saw in the case of com- 
plex u. Note that for the special case of ^ = 0, the exact 
solution (cf. (p2)) ) is proportional to 5(cos0). The scale 
parameter of the Gauchy distribution will be 



7' 



gr,/-D J. 

WD)r- 



(34) 



Only the transients of the distribution add to the phase of 
6, and numerical simulation shows that they are enough 
to again make it uniform. 



IV. DISCUSSION 

In the context of the NLSE, the scattering data of the 
ZS operator uniquely determine the solitonic excitations 
we get in the optical fiber if we feed one end with a delta- 
correlated Gaussian signal. Even though the informa- 
tional contents of the Gaussian signal and its solitonic 
spectrum arc the same, it is easier to consider the effect 
of amplifier noise in the domain of the scattering data. 
For example, a small amount of amplifier noise will ran- 
domly shift each eigenvalue z by a small amount, while 
making large changes in the output signal [H| ■ The effect 
of this noise is important to analyze, in order to calcu- 
late the ultimate information capacity limits through op- 
tical fibers. In principle, to find the capacity one needs 
to optimize over input signal distributions, which is a 
formidable task. Instead, in this paper we start with a 
given input distribution and calculate the corresponding 
density of states and the corresponding distribution of 
scattering data b^ ■ We leave the analysis of the effects of 
noise on the spectrum for a future publication. 
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APPENDIX A: THOULESS FORMULA 

The proof of the Thouless formula for the ZS eigen- 
problcm proceeds similarly to the proofs in [2^, [2§| . We 
consider the system of equations 1^ on the interval [0, T] . 
Let *'(i), *'■(*) be two independent solutions of ^ that 
satisfy the conditions: 



*'fO) = ( ^^^'^ 



^ ^'>-[r2iT) ) 



(Al) 
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We will need to combine this pair of initial and final 
conditions into a set of boundary conditions for the 
eigenstates, and for this we let each of them be a one- 
parameter family of initial (final) conditions to avoid over- 
determining the problem. This means that '012(0) are 
not chosen independently, but satisfy a single linear rela- 
tion. The same goes for ■(/i[2(2^)- The Wronskian of the 
two solutions, W = V'i'02 ~ V'SV'i is constant. Taking the 
derivative of ^ with respect to z, we obtain an equation 
for dz'^^{t) whose solution can be written in terms of a 
matrix Green function G': 



d_ 
dz 



dt'G\t,t')<b\t') (A2) 

^0 



with 



G\t,t') 



Here cri 



t > t' (A3) 
i < t' 



1 



1 

isfies the initial conditions 



. The matrix Green function G sat- 



G\0,t') = 0, §^G\t,t') 



0. 



(A4) 



4g|. 
Ml 



where 'Texp' denotes the path-ordered exponential 
Because the trace of V vanishes, we have det S = 1 
Therefore, if we denote the two eigenvalues of S{t, 0) by 
s±{t), with |s+(t)| > \s-{t)\ = |s+(t)|-\ we have 



lim - In I s+(T) I, 



|._(r)^_(T)r 

(AlO) 



independent of the initial condition, given that the coef- 
ficient of the exponentially increasing solution does not 
vanish, which is the case with probability one in the limit 
of large T. 

Now, since ^'^(T) depends only on z* and not on z, 



taking the derivative of (|A6p with respect to the latter 
we get 



—— = lim — - 



(All) 



The quantity inside the average can be computed from 
[1: 



t=o 



We also define another Green function. 



G(t,0 



r*'-(t) (<Tl*'(t'))^, t > t' 



(A5) 



which satisfies the conditions (|A1[) (taken together as 
boundary conditions) and will determine the density of 
states. The Lyapunov exponent can be expressed as(48[: 



(A6) 



Before going any further, we must note that the value 
of K is, with probability one, independent of the initial 
conditions satisfied by 4*' (the argument is very similar 
to that for the FGM [13]). To see this, we rewrite the 
system of equations ^ as: 



d 



1p2 



V 



"02 



z —u 
—u —z 



(A7) 



*'^(r)a3*'(T) 
*'^(T)*'(r) 



.*'^(T)*'-(T) j^dt' {2ijA{t'W^{t')) 



w 



[ dt'TTG{t',t') (A12) 
Jo 



The first term on the right hand side of (|A12|) is almost 
surely 0(1) in the limit of large T and so does not con- 
tribute to the average. We are thus left with 



|.^lim2,/£..'TrG(t',0) (A13) 



Taking the derivative of (|A13p with respect to z* and 
using the relationfl^l p{z,z*) ~ (TrG) wc arrive 

at the Thouless formula for the density of states: 



We can formally write the solution for 4^': 
^\t) = S'(t,0)*'(0) 
5(^,0) =rexp 



-I I dt'V{t') 
'o 



(A8) 
(A9) 



p[z,z ) = - 



ndz*dz 2Tr [oe ^ dr^^ ' ^ ' 
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Note that in the analysis of light propagation through 
optical fibers the roles of space and time have been 
exchanged, as compared to the traditional nonlinear 
Schrodinger equation: Propagation down the fiber has 
taken the place of time, while the traditional space vari- 
able has been replaced by the retarded time measured in 
the frame moving along the fiber with the group velocity. 
In the case of a single localized eigenstate with z — ^ + ir] 
and b = |&|e"'', the corresponding soliton has amplitude ^, 



— i^lp and initial phase 



velocity 2»7, initial "position" to 

00 = - TT 01 . 

Regarding the initial positions (in time) of the solitons, 
this means that the solitons will be uniformly distributed, 
independently of velocity or amplitude in the limit T 
00. This is expected, given the translation symmetry of 
the initial condition. 

The average over the Gaussian enesmble is included here 
so that our final formula refers to the average DOS. It is 
not needed for the definition of k. 

Note that V is Hermitian (and S{t,0) is unitary) when 
z is real, and so the Lyapunov exponent is always zero 
in that case, irrespective of u. Therefore, any eigenstates 
on the real axis are extended, corresponding to scattering 
states of the ZS eigenproblem. 



